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Abstract 

We present the complete electroweak one-loop corrections to the partial widths for two-body 
decays of a chargino (neutralino) into a VT-boson and a neutralino (chargino). We perform 
the calculation for the minimal and the next-to-minimal supersymmetric standard model 
using an on-shell renormalization scheme. Particular attention is paid to the question of 
gauge invariance which is achieved using the so-called pinch technique. Furthermore we 
show that these corrections show a strong parameter dependence and usually are in the 
range of 1-10 percent if the neutralino involved is a higgsino or wino like state. However, in 
case of a bino-like or a singlino-like neutralino the corrections can go up to 50% and more. 
Moreover we present the public program CNNDecays performing these calculations. 
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1. Introduction 

A new energy frontier has been opened with the start of the LHC, namely the exploration 
of the tera-scale. An important part of the experimental program of the LHC is the search 
for supersymmetric particles (J, 0, S]- ^ they are indeed found, the determination of the 
underlying parameters will then be among the main tasks. Here the properties of neutralinos 
and charginos will play an important role as they occur in various steps in the cascade decays 
of squarks and gluinos which can be copiously produced at the LHC. These investigations 
will then most likely be complemented and continued at a prospective future linear collider 
(ILC) 0, El and/or at a multi-TeV collider such as CLIC 

The two-body decay of a chargino (neutralino) into a VT-boson and a neutralino (chargino) 
is one of the most important modes if kinematically allowed js), Q, The determination 
of the corresponding branching ratio can give important information on the nature of the 
chargino and neutralino involved and therefore one will need at least an ILC for its precise 
value. This is one of the reasons to investigate the one-loop corrections to this decay mode. 
We will see below that the corrections can be as large as the LO decay width in certain 
regions of parameters space. A second important aspect is the question how to perform the 
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complete electroweak corrections in a gauge invariant way. These decay modes offer an ideal 
play ground to tackle the technical questions involved. 

We will use an on-shell renormalization scheme for masses and couplings. The on-shell 
renormalization of the mass matrices of charginos and neutralinos in the minimal super- 
symmetric standard model (MSSM) was first discussed in 0, [h| and extensively addressed 



in 
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Moreover, electroweak corrections for these decays taking only into account third 
generation squarks and quarks within the MSSM have been calculated in [12| . Corrections 
of up to 10% have been found there. Note, that for these corrections gauge invariance is not 
an issue. 

We will work in the general linear i?^-gauge. An important aspect will be the question 
how to renormalize the mixing matrices of charginos and neutralinos, since it is already 
known from the Standard Model that the on-shell renormalization prescription according 



to 



13l . Il4| leads to a gauge dependent CKM-matrix [15( which in turn implies a gauge 



dependence of processes like t — > Wb at next-to-leading order level [16 



Since then, some papers proposed different solutions to the problem addressed above: 
whereas [l7| argued that missing absorptive parts due to the unstable nature of the external 
particles have to be included in the calculation, [l6[ proposed a method on how to construct 



a gauge invariant counterterm for the mixing matrices inspired by the pinch technique |18 
which defines gauge independent form factors for gauge bosons. Another perspective is 



presented in [19|, where the gauge- variant on-shell renormalized mixing matrix is related to 
a gauge independent one in a generalized MS scheme of renormalization. 

For the special case of the CKM matrix different methods were discussed in the literature: 



15| proposed to set the momenta of the non-diagonal entries of the quark self-energies to 



zero. 



20l . l2ll |22| suggested new variants of renormalization in particular for the mixing 



matrices themselves partially based on physical processes which allows one to get gauge 



23l . [24| proposed 



independent decay widths. For lepton or neutrino mass matrices also 
useful renormalization conditions allowing gauge independent results. 

In this paper we will apply the technique presented in [l6j]. The technical features are 
the same for the MSSM and the next-to-minimal one (NMSSM). As for the processes under 
consideration the MSSM is a trivial sub-class of the NMSSM in the limit of heavy singlet 
states we will perform the calculation right from the start in the NMSSM. The paper is 
organized as follows: in Section [2] we fix our notation for both models and in Section [3] we 
summarize the tree-level results for the mass matrices of charginos and neutralinos as well 
as the tree-level decay widths for \t ~~ >* Xk^ ± an d Xi ~^ xfW T . The virtual one-loop 
contributions and real corrections including the question of gauge invariance are discussed 
in Section |H We present our numerical result in Section and in Section O our conclusions. 



Appendix A contains the generic formulas for the matrix element contributions of the 



vertex corrections. In |Appendix B| we give the formulas for the the real corrections which 
are non-factorizable as the higgsinos form a vector-like SU(2) gauge state. The usage of the 
i?g-gauge implies that all possible derivatives of the two-point one-loop functions show up in 
the calculation. Some of them are to our knowledge not available in analytical form in the 



literature and, thus, we provide them for completeness in Appendix C| In |Appendix D 



we present the program CNNDecays for the numerical evaluation of these corrections. 
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2. The models 



In this section we fix the notation for the models considered, the MSSM 251 and the 



NMSSM [26|]. Detailed formulas for mass-matrices and couplings can be found in [25l. 127 



Both model have in common the Yukawa part of the matter superfields coupled to the 
ST/(2)-doublet Higgs fields 

W y = e ab (YtiWfil + YSQlDjHS + Y?L%Hl 



where e ab is the complete antisymmetric SU(2) tensor with eu = 1. In case of the MSSM 
the /^-term is added and the total superpotential reads as 

Wmssm = Wy- ^H d H u , (2) 
whereas in the NMSSM the Higgs doublets are coupled to a gauge singlet superfield S 

Wnmssm = Wy- XSH d H u + IkSSS . (3) 



Wy - \SH d H u + -kSSS 

If the scalar component of the gauge singlet S gets a vacuum expectation value 4=u 5 an 
effective //-term is generated 



(4) 



In this way one gets fi naturally of the order of the electroweak scale [26|. To complete 
the models one has to add the soft-SUSY breaking terms, where the part common to both 
models reads aa3: 



V} oft =mfQTQ) + ">rr,r; + mp,^ + mflfL) + mfE.E* 



i2~ ~ 



2 ~ ~ 



;2 ~ 



;2~ ~. 



1 

+ 2 



M^B" + M 2 W C W C + M 3 g d g d + h.c. 



+ tab 



rjQ«U*H b u + T%Q\D*H a d + rfL\E*H a d + h.c. 



(5) 



with a,b = 1,2 as summation indices. The complete soft-SUSY breaking potential of the 
MSSM is then given by 

V™f? M = Vl ft -B^ ab H a d H b u (6) 

and the NMSSM one by 



taNMSSM 
v soft 



V s \ ft + m 2 s SS* - [e ab T x SH a d H b u + h.c] + 



-T K SSS + h.c. 



(7) 



1 We use the notation of [H, E 



3. Tree-level results 



Here we collect the tree-level results which form the basis for the one-loop corrections 
considered in the subsequent sections. 

3.1. Masses of Charginos and Neutralinos 

Using Weyl spinors in the basis (ip~) T = (W~,H~[), i^ + ) T = (W + ,H+) the chargino 
mass matrix is given in both models by 



M c 



M 2 j= 2 gv u 



T2 gv d [i 

The mass eigenstates are obtained by two unitary matrices U and V 



F+ = V lt ^t and F: 
such that the mass eigenvalues are given by: 

Diag(m~+,m~+) = 



U*M r V 



-1 



(9) 



(10) 



m-+ < m~+ 

Xl A2 



and xt is the Dirac spinor given by 



X.i 



FI 



(11) 



For the neutralinos we obtain in the basis^l (?p ) T = (B°, W®, H^, H®, S) the mass matrix 



( 


Mi 





-\g'v d 


\g'v u 





\ 







M 2 


\gv d 


-\gv u 









-\g'v d 


\gvd 







"75 A ^ 






\g'v u 


-\gv u 


-// 









\ 








"73 A ^ 


"73 A ^ 


\/2kvs 


J 



(12) 



where in the NMSSM /i is given by (J4]). In case of the MSSM one has to omit the last 
column and the last row in Equation ( TT2|) . The mass eigenstates are obtained via 



F? = 



(13) 



from the gauge eigenstates if?®, where the unitary matrix M diagonalizes the mass matrix 
M. n as 



Diag (m^o, 



Af*M„Af j 



(14) 



2 S is of course only present in the NMSSM. 
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where j = 4(5) in case of the MSSM (NMSSM). The masses are ordered as m^o < m^o for 
i < j and we have chosen Af such that all mass eigenvalues are positive. This implies that 
in general the matrix M is complex. Useful checks can be performed if one goes in case of 
real parameters to a basis where M is real. Then one or more of the mass eigenvalues get 
negative and the mass ordering applies to their moduli. For completeness we note that Xi 
is the 4-component spinor given by 

X? = ( S ) • (15) 



3.2. Decay widths 

The partial widths for the decays x°j ~> xfW T and xf —> X°jW T are obtained from the 
following interaction Lagrangian: 

C = X?7 M (P L Llj + PrO RIj ) x-W- + h.c. . (16) 
In both models the couplings are given by: 

Olij = - gU; 2 U a - O m = -gVZATp + ^/Wfo ■ (17) 

The widths have the form 

^ = 16^M<<)^El^| 2 , (18) 

* pol 

where m 8 (m ) is the mass of the mother (daughter) particle and M T is the tree-level matrix 
element. Explicitly they are given by 



with 



x((|C^| 2 + \0 Rl3 \ 2 )f(m 2 ,m 2 ,m 2 w )-6Re(0 Llj O Rlj )m,m^ (19) 
x((|0 Llfe | 2 + \0 Rlk \ 2 )f(m 2 ,ml,m 2 w )-me(0 Llk O Rlk )m l m k ^ (20) 



f(x,y,z) = I(x + y )-* + ^j£ (21) 
k(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz . (22) 
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F S V 




Figure 1: Generic self-energy diagrams for gauge bosons. 



S V 




F F 



Figure 2: Generic self-energy diagrams for fermions. 



4. One loop results 

We discuss the one-loop corrections in the context of the NMSSM as it contains the 
MSSM in the limit of very heavy gauge singlet states. There are no technical differences for 
the decays considered as the gauge parts of both models coincides. The one-loop corrected 
(renormalized) amplitude(s) M can be expressed as 

M = M T + AM = M T + M v + M wv , (23) 

where My contains the vertex corrections and Myyy is the sum of the wave-function cor- 
rections and counterterms. Moreover we consider finite mass shifts m~±o — > ml± for the 

charginos and neutralinos at the one-loop level, since the number of physical parameters is 
smaller than the number of masses. 



4-1- Wave-function and coupling renormalization 

General formulas for the wave-function renormalization of supersymmetric particles can 



be found e.g. in j30[. Here we summarize the most important points related to the decay 
under consideration and discuss in particular the problem of gauge invariance. We work in 
a linear gauge given by 



C 



■^7 



_L(i^)2 _ J_ F + F - 
2£z <iw 



(24) 
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with 



F ± = d"W^ q= iM w i w G ± (25) 



F z = 0% - M Z £ Z G (26) 
F 7 = d^A^ . (27) 

The inverse propagator of the vector bosons is given at tree-level by 

-ig^{k 2 -ml) + i(l- ±^ = -ig^{k 2 - m\) - ig? (m 2 v - ^ (28) 
with V = W, Z and the transverse and longitudinal projectors are defined as: 

& = tr - & = ^- (29) 

The multiplicative renormalization of the bare parameters is done in the following form: 

m v0 — > Z m m v = (1 + 5Z m )m v = m v + 5m v with 5m v = m v 5Z m 
V° -> ZyVp = (1 + 15Z V )V„ 

ZVO -+ Z-^y = (1 + 5Z iv Y^y (30) 

The Z^s are fixed by the conditions that propagators mixing the Goldstone bosons with the 
vector bosons vanish [3l[ and will not be discussed further here as they are not needed in 
the subsequent discussion. The vector boson propagator then reads at the one-loop level 

T^(A; 2 ) = _ <{k 2 _ m 2 y) + if u _ £) _ ig ^t v T {k 2 ) - y?±l{#) (31) 



with 



t^(k 2 ) = T^{k 2 ) + k 2 5Zy-m 2 v ~5m 2 y 



Yy L (k 2 ) = Z v L (k 2 ) + tm v 8Zv + &m v -k 2 (5Zi: + 5Zv) , (32) 

where T^f(k 2 ) and E^(fc 2 ) are the contributions obtained from the one-loop graphs whose 
generic structures are shown in Figure [TJ 

We require an on-shell renormalization for the vector bosons, e.g. that the pole of the 
propagator occurs at the physical mass and that its residuum is one. This yields the condi- 
tions 

Return 2 ,) = (33) 



~*v f 2 . <9Re£^(fc 2 ) 
ReS' T (m 2 / ) := 

8 



. (34) 



k 2 =m v 



Note, that the Re applies only to the one-loop functions but not to the couplings. Finally 
one gets: 

5rriy = ReSy(my), 5Z V = — ReE'^(my) (35) 

Moreover one has to fix the renormalization of the couplings e and g and the cos of the 
Weinberg angle cos 9w'- 

e 5Z e e = (1 + 5Z e )e = e + 5e (36) 
g -> 5Z g g = (1 + 5Z g )g = g + 5g (37) 
cos 9w — > cos 9w + 5 cos 6w (38) 

Not all of these quantities are independent of the mass renormalization and one obtains [l 

. . 1 .. f SmL 5m 2 7 \ . 1 . . 

o cos 6V = - cos v w I — 2 2~ == ' > isin% = - - — — dcos9 w (39) 

z y m w m z J tan u\y 

« e= W)_^E^(0) =► «,= (f2.-^t), (40) 

2 777,^ \ Sill u\y J 

In case of charginos and neutralinos the mixing effects have to be taken into account. 
The fact, that both sectors together have less parameters than mass eigenstates, leads to 
additional finite shifts and has been discussed in detail for the MSSM in @, Ed, Q being 
relevant for the calculation of NLO masses. Before addressing this topic we first present the 
principal idea of the wave-function renormalization for mixed fermions and discuss gauge 
invariance in all detail. For the wave-function and mass renormalization one needs 

X? -> (<% + \5ZfP L + \5ZfP R ) x°j (41) 
Xt (<% + \5ZfP L + \5Z± R P R ) xf (42) 
rrifi — >■ rrifi + 5m fi . (43) 

It is well known that the wave-function renormalization constants of the neutralinos and 
charginos satisfy the following relations 

5Zf = 5Zf\ 5Z± L = 5Z~ R *, 5Z~ 3 L = 5Z± R * . (44) 



As already pointed out in [13J the combination of tree-level mixing with non-diagonal field 
renormalization gives anti-hermitian parts, which however can be canceled by introducing 
counterterms for the mixing matrices themselves: 

Mij -> Afij + 5Mij with SAfij = i (5Z°£ - 5Z° k R ) Af kj (45) 

U« -+ Uy + SUy with 811* = ~ (5Z+ R * - 5Z+ R ) U kj (46) 
V ..^ V .. + 5V .. with 5V .. = \ (sz+l _ 6z +l*) Vk . (47) 



Here we recall the basic steps to renormalize a system of Dirac fermions. The case of 
Majorana fermions can be obtained from this by noting that the Majorana nature induces 
relationships between Z-j and Z R . As in case of the vector bosons one starts with the inverse 
propagator given by 



if£ (p) = ifiyO* - m fi ) + i [j> (PlH- (p 2 ) + P R £i R (p 2 )) 

+p L % SL (p 2 ) + p R % SR (p 2 ) 



with 



% L (P 2 ) 


= n-{p 2 ) + 


% R (P 2 ) 




% SL (P 2 ) 


= sr (p 2 ) - 


n? R {p 2 ) 





(48) 

(49) 
(50) 
(51) 
(52) 



The corresponding generic one-loop contributions are shown in Figure [2j The on-shell con- 
ditions read as 



Ui(p)Ret{Ap) 



ReK^ujip) 



1 9 









__ \fT"A// $ + m fi 

hm n u^ReT^p) 

P " L fi 



p 2 —>m 2 „ . 



Uiip) 



Si 



2 Um 2 T1 2 +7 I / 2 Re %(P) U i(P) = U i(P) 
P 2 ^rn 2 ft p —m fi 



(53) 
(54) 



which in turn lead to the counterterm for the masses rrifi 
1 



5m 



fi 



mfiRBEt (m fi ) + m/;ReS£ (m fi ) + ReE^ (m fi ) + ReE 



(55) 



and the following wave-function renormalization constants 



fL 



5Z 



8Zl R 



-Re 



-Re 



E^m^ + m^fE'^KJ + E 



+m /4 (E'f i K i ) + E 



+m /i (E'ir(ml : )+S'ir(m 



fit 



"fi. 



(56) 



(57) 
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for the diagonal entries, whereas for the off-diagonal ones one finds: 



m% - m 2 fj 



m% - m% 



+m/ i ReSf J - L (mj j ) + m fj ReE 



fSR. 



m 



mfiTnfjReTi 



+mfj~ReT, 



l SL ( m }j. 



+ m/jReS 



fSR. 



m 



Si' 



(5? 



(59) 



4-2. Gauge invariance 

Before discussing the finite mass shifts for neutralinos and charginos and the question 
how to construct Myyy in the next subsections we have to clarify a subtle point with re- 
spect to gauge invariance. The counterterms given in Equations f l56|) - fl59l) cancel the UV 
divergences and avoid anti-hermitian parts in the Lagrangian. However, this does not nec- 
essarily imply gauge invariance for the process as well as the one-loop masses. To achieve 
this we use the method proposed in 16[ which is inspired by the pinch technique 18j, which 
defines gauge independent form factors for gauge bosons. We are aware that it has one 
weak point, namely a dependence on the choice of the gauge fixing for the mixing matrix 
counterterm. Its advantage is that it is model independent and does not rely on the details 
of the renormalization of physical parameters. 

The basic idea of this method is to treat the gauge dependence of the counterterms for the 
wave-function renormalization constants differently from the ones for the mixing matrices. 
Therefore we calculate two variants of wave-function renormalization constants, namely 
in case of the neutralinos SZ^ L , SZff 1 for arbitrary values of £y and SZ^ L , SZ^f 1 for £v = 1 



('t Hooft-Feynman gauge). The same is done for the wave-function renormalization constants 
of the charginos which we denote by SZ^ L , 8Z^ and SZ^ L , 5Z^ R . The counterterms for 
the mixing matrices are calculated via the wave-function renormalization constants in the 
't Hooft-Feynman gauge 



5Z™ ) M 



SUa = \ (<^T - ^) U kj , SVa = \ (5Zf k L - 62$*) V kj 



(60) 
(61) 



whereas for the wave-function renormalization constants as given in general form in Equa- 
tions (I5T)]) - (159"|) we use the full gauge dependent form. This splitting forces one to include 
the additional contributions from the tadpole graphs with the Goldstone bosons as shown 
in Figure [3] A few remarks are in order here: for £y = 1 the tadpole graph contributions to 
the wave-function renormalization exactly cancel their contributions to 5J\f, 5U and 5V in 
case of the considered processes. For (7 ^ 1 they cancel also £y-dependent parts of other 
one-loop contributions. Moreover, they only contribute to the ultraviolet divergence for 
^ 1 but of course the rearrangement is such that the total ultraviolet divergence cancels. 
The infrared divergence in case of the processes is not affected at all by this mechanism, 
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/ / 

It , Lx I At , Lx I 

/ / 



I I I I 



4- 



X/ Xj X/ X,j 

Figure 3: Tadpole contributions induced by Goldstone bosons G , G° ensuring gauge invariance. 



since the photon only contributes to the diagonal entries of fermionic self-energies. Last 
but not least we note that all these statements hold as well for the calculation of gauge 
invariant, ultraviolet and infrared finite masses ml± in the next section, implying that we 
use gauge invariant counterterms for the mixing matrices including the described Goldstone 
contributions. 

4-3. On-shell masses of Neutralinos and Charginos 

A special feature of the chargino/neutralino sector is, that the number of parameters 
is lower than the number of imposed on-shell conditions. This implies finite corrections to 
the tree-level masses of neutralinos and charginos. We closely follow [9] and start with the 
one- loop contributions to neutralino masses 5A4.® and chargino masses 5Mf 

{^■^n)ij - $ {N T M® dia N).. (62) 

n,l 

mt) ij = HU T MZ lia y) i . (63) 
= [ 5u m (Ml dia ) nl V v + U ni {5M% dia ) nl V v + U ni {M® dm )jV lj _ 

n,l 

with the diagonalized mass matrices Ai dia and their counterterms SAi dia : 

(Ml dm ) nl = 5 nl m^, (SMl dm ) nl = S nl Sm^ (64) 
(M® dia ) nl = 6 nl m^, (5M® dia ) nl = 5 nl 5m^ (65) 

The bare mass matrices of the neutralinos and charginos can now be expressed as the correct 
on-shell mass matrix with the corrections ( 162]) and (163]) or via the tree-level mass matrix 
expressed in physical parameters together with the renormalization constants of those: 

AC = Ml c + 8Ml c = M n , c + 5M n , c (66) 



12 



Therefore the relations between tree-level and one-loop mass matrices take the form: 

AC, C = M n , c + SM n , c - 6M% e =: M n , c + AM n , c (67) 

In the following we will define the model-dependent physical parameters, write down the 
renormalization of the mass matrices 5Ai n ,c and identify the renormalization constants of 
the physical parameters, which should be fixed in the neutralino or chargino sector. Some 
physical parameters, namely dmw, 5m z and thus 6 cos 9w are fixed in the gauge boson sector 
and S tan 8 in the Higgs sector. In particular for 5 tan B we take the DR renormalization 
fjil ]. such that UV divergences in the masses and the considered process cancel 

5 tan 8 1 



r A (3Tr(y d y d t) - 3Ti(Y u Yj) + Tr(Ytf)) 



(68) 



tan 8 32vr 2 

with A = - — 7 + ln(47r). Note, that this choice maintains also the gauge invariance of 
masses and the considered decay widths. 

4.3.1. MSSM 

In case of the MSSM we follow [9]. The tree- level neutralino mass matrix in (Fl2|) can be 
written as 

/ Mi — m z sin 9 W cos (3 m z sin W sin [3 \ 

M 2 mz cos 9w cos (3 — mz cos 9w sin (3 

-n 

\ sym. J 

and the chargino mass matrix in (jHJ) as 

M 2 \/2mw sin (3^ 

^\/2m w cos(3 /j, 

The variation of all given entries of the tree-level mass matrix leads to 

5Mx 



M n 



(69) 



(70) 



5M1 1 



5M\f 
SMI 4 
5Mf 



SMI" 



5M X 



-M 11 

Mi n 



—5(mz sin 9w cos (3) 

S(mz sin 9w sin (3) = 
SM 2 



Smz 5 sin 9w 5 cos 8 
+ . -„ + 



m z 



Smz 
m z 



+ 



sin^ 
5 sin 9w 



sin 6 1 



+ 



w 



cos /3 
5 sin B 
sin /3 



M^ 3 
-Mi 4 



SMo 



-M 



22 



5M 



31 



M 2 

5(m^ cos 6*iy cos (3) = 

—5{mz cos 6*w sin 3) 
5fi 



5mz 

m z 
Sm z 
m z 



S cos 9w 8 cos f3 

cos 9w 

S cos 9w 
.{ y JL 

cos 9w 



+ 



cos 6 
S sin B 
sin/3 



M 23 



M 



2 I 



-5fi 



-M 



3d 



(71) 
(72) 
(73) 
(74) 
(75) 
(76) 
(77) 
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whereas all the other variations 5M 12 = SM 21 = SM^ 
The corrections in the chargino mass matrix read 

5M 2 



5M 



14 



6M? 



5M? 
SMf 
5Mf 



-M 22 



y/2S(mw sin (3) 
\/25(m w cos/3) 



5m w S sin (3 
m w sin (3 

5mw 5 cos (3 
mw cos (3 



M 
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A* 



necessarily vanish. 

(78) 
(79) 
(80) 
(81) 



21 



We will fix M 2 and \x in the chargino sector, whereas M\ is fixed in the neutralino sector by 
imposing the conditions 



AM 



ii 



AM 



22 



AM 



ii 







resulting in: 



8M l = SM 



©ii 



5M 2 = 5M 



©ii 



8n = 5 Ml 



;*;22 



(82) 



33) 



For all the remaining entries of the neutralino and chargino mass matrices finite shifts AM n ,c 
have to be taken into account. 



4.3.2. NMSSM 

Defining the additional angles and parameters 

tan(3 s = — , n = —=\v s , 
v u V2 

the tree-level mass matrix of the neutralinos is given by 
/M 1 



m s 



\/2kvs 



54) 



M n 



— m z sin 9 W cos (3 mz sin 9 W sin j3 

M 2 mz cos 9 W cos (3 —m z cos 9w sin (3 

-fi 
sym. 








V 



tan /3s 
tan tan /9g 

m 5 



i5) 



and the chargino mass matrix is equal to the one in the MSSM. Beside the variations (j71[) - 
(|77p and (JT8"]) - (j8~I]) already present in the MSSM one has in addition 



SMf 



5M 



55 



tan (3s 

tan (3 tan (3s 
8m s 



5(jl 5 tan (3s 



M 



35 



/i tan (3s 

5fi 5 tan (3 5 tan (3s 
fi tan /3 tan (3s 



M 



45 



(86) 
(87) 



55 



m 5 
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whereas all the other variations 5M 12 = 5M 15 = 5M 2 } = SMf = SMf = SM 44 = 
necessarily vanish. Similar to the MSSM we will fix M2 and /i in the chargino sector, 
whereas Mi, tan (3$ and rris are fixed in the neutralino sector by imposing the following 
conditions 

AMI 1 = AMf = A ^ X n = AMf = AM 5 n 5 = » ( 89 ) 

which results in: 

8M1 = 5M® 11 , 5M 2 = 5Mf\ 5(i = 5Mf 2 (90) 
tan 2 f} s ( \ , 35 1 34 \ — @55 



5t an/3 s = ^ 5Mf 5 - —-5M®* 4 , 5m 5 = (91) 

/i \ tanp 5 J 

For all the remaining entries of the neutralino and chargino mass matrices finite shifts AM n ,c 
have to be taken into account. Note, that one could also fix £ tan /3s in the Higgs sector. 

4-3.3. Effect on the considered processes 

With the procedure introduced in the last sections we can calculate one-loop on-shell 
neutralino and chargino masses for the MSSM and the NMSSM, namely by combining 
the full one loop corrections SM® c with the counterterms SM n ,c obtained according to 
Equation ( 1ST)) . This results in the one- loop mass matrix M® c , which diagonalizations lead 
to one-loop neutralino m l -o and chargino masses mi£ and mixing matrices at the one- 

loop level M lL , U 1L , V 1L . Note, that these masses are UV and IR finite as well as gauge 
independent, if one takes into account the gauge independent renormalization of the mixing 
matrices f )62|) and f )63|) . 

Instead of having the diagonal counterterm for the masses in f )4"3")) we have to replace 

Sijmfio ->• Sijmfi + SMijP L + SM^Pr , (92) 

where 5M = D* R 5MD^ L with 5M being the physical mass counterterm and Dl,Dr being 
the rotation matrices, which diagonalize the tree-level mass matrix Mdia. = D r MDl in 



the notation of 11]. With this counterterm one gets contributions to the non-diagonal wave 



function renormalization constants: 



5Z[ L 



13 m % - m % 



m 



Re£{/(my + m /i m /j ReEf/(mJ i ) 



fj 



+m /i Re£{/ i (m 2 v ) + m fj ReZ{f R {m 2 fj ) - m i 8M ij - rrijSM^ (93) 



szf R 



2 



13 m % ~ m )j 



m /i m /i ReE^(mJ j ) + m 2 fj ReE{ R (m 2 fj ) 



+m fj ReE{f L (m 2 f] ) + m fi ReZ{f R (m 2 fj ) - mj 5M i3 - rmSM^ (94) 

However it turns out, that these additional contributions are cancelled by the contributions 
to 5J\f, SU and 8V, which also have to be calculated using the new wave-function renormal- 
ization constants. This implies that the reduced number of physical parameters only has an 
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impact on the calculated neutralino and chargino masses, but not directly on the considered 
processes itself. 

The fact, that one has finite shifts to the on-shell masses at the one- loop level leads to a 
complication: to obtain infrared finite results one has to take into account the emission of an 
additional photon as discussed in Section 14.61 Using one-loop corrected masses in these pro- 
cesses requires that one also has to use one-loop corrected masses in the virtual corrections. 
Formally this results in differences which are of higher order and indeed numerically these 
differences are rather small. However, the use of one-loop masses in the virtual corrections 
leads to a spurious gauge dependence which is of two-loop order and which eventually gets 
canceled by taking into account two-loop corrections. We have checked numerically in the 
examples presented in Section that indeed this residual gauge dependence is very small. 

4-4- Counterterm for the considered processes 

Plugging everything together we get from the tree-level interaction in Equation (Tl6|) the 
counterterm Lagrangian 



5C D xiY Pl 



1 1 5 1 2 

5 °Lij + -0 Lij SZ W + - L ikSZ° Lkj + - 5Z Lki°Lkj 



k=l 



k=l 



(95) 



+Pr 



1 1 1 

SO Rij + -0 Rij 5Z w + - °Rik 5Z °Rkj + o X] 5Z Rki°Rkj 



k=l 



k=l 



with 



SO Lij = - (5gAf* 2 U a + g6Af* 2 U a + gN* 2 5U a ) 

- ±= {5gN* 3 U i2 + g5M* 3 U i2 + gAf* 3 5U i2 ) 
SO mj = - (6gV£Af j2 + g5V*Af j2 + gV£6N j2 ) 



+ -j= (6gV£Af j4 + gSV^fj, + gV^Af^) 



(96) 



(97) 



Therefore we define 



/ 1 i 5 i 2 \ 

5A Lij = i 60 LlJ + -0 LlJ 5Z w + Ou k SZ° Lkj + - SZ L * kl O Lkj 

V k=l k=l J 

SA mj = i 50 Rij + -0 Rij 5Z w + 2 °R^ Z Rkj + ^ Yl 6Z Rki°Rki I 

\ ~ k=l k=l J 



(98) 
(99) 



and obtain 



M wv = u( Pl )r (P L 5A Llj + P R 5A Rlj ) u{k)el{p 2 ) 



(100) 
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d) 




F 



F„ 



Vo 



Figure 4: Generic Vertex NLO corrections. 



4-5. Vertex corrections 

The next piece for the one-loop corrections are the vertex corrections. In Figure H] we 
depict the 6 generic contributions to My. In the Feynman diagrams a) and b) fermions and 
sfermions contribute as well as charginos/neutralinos together with the neutral and charged 
Higgs bosons. In diagrams c) and f) only charginos, neutralinos and the vector bosons 
(including the photon) contribute whereas in diagrams d) and e) there are in addition the 
charged and neutral Higgs bosons as well as the Goldstone bosons. The individual contri- 



butions from the diagrams in Figure H] to the matrix element My are given in Appendix A 
for the 't Hooft-Feynman gauge £y = 1. The general case ^ 1 leads to rather lengthy 
formulas which are included in the program CNNDecays ji^. Putting everything together 
we obtain for the width at the one-loop order including virtual corrections only 



r lL = pO + ^W,<<) 1 y 2Re ((MwM *) + ( My M* )) . (101) 
lo7rmf 2 ^— ' 

pol 

This result is UV finite, gauge independent and also not dependent on the renormalization 
scale Q, which cancels between the contributions from My and Mwv- However we still have 
to address the IR infiniteness in the next section. 

4-6. Real corrections 

Last but not least one has to take into account the real corrections to cancel the infrared 
divergences occurring due the photon contributions in the self-energy diagrams in Figures [1] 
and [2] or the vertex corrections in Figure HI Technically this is achieved by introducing 
a small photon mass for both, the virtual corrections and the real photon emission. The 
corresponding graphs are shown in Figures and O We calculate the real bremsstrahlung 



corrections according to [14| and collect the rather lengthy expression in Appendix B We 
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Figure 5: Real corrections for Xj ~^ X; F W /± 7- 




Figure 6: Real corrections for xf — > x?W /± 7 



note, that (i) due to the presence of left and right couplings, see (Fl6|) . these corrections do 
not factorize. (ii) The gauge dependence of the graph with the Goldstone boson G + cancels 
the one with the virtual M^-boson implying that the real corrections are gauge independent. 

Denoting the width due to real bremsstrahlung by T R we obtain for the IR-finite decay 
width at next-to-leading order 



5. Numerical results 

Let us discuss a general feature before presenting our numerical results. Inspecting 
the couplings given in f[T7|) one sees that only the wino and higgsino components of the 
neutralinos couple but not the bino or the singlino one. This implies a small tree-level width 
if the neutralino involved is either mainly bino or singlino which potentially results in large 
corrections to the widths considered as we will see below. 

In the subsequent section we first present our results for the NMSSM benchmark points 



given in [35|, |36[ . Afterwards we consider cases where the neutralino is nearly a pure bino or 



singlino. There we demonstrate the additional corrections not considered so far as important 



as the squark/quark contributions considered in the literature 12 



5.1. NLO masses and decay widths for benchmark scenarios 

In Table [T] we recall the most important parameters of the benchmark scenarios used in 
our analysis based on 35| and j36j. They are minimal supergravity (mSUGRA) [351] and 
gauge mediated SUSY breaking (GMSB) [36[ scenarios. The parameter points can also be 
found within NMSSM-Tools El. 
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moUGKA 1 


mbUGKA 6 


nTT/lT) A A 


GUI Mq (in GeV) 
scale M1/2 (in GeV) 
A a (in GeV) 


180 
500 
-1500 


178 
500 
-1500 


TO A 

780 
775 

-2250 


SUSY tan/3(m z ) 
scale fj, (in GeV) 

A, K 

A X ,A K (in GeV) 


10 
969 
0.10,0.10 
-959,-1.6 


10 

938 
0.40,0.30 
-616,-11 


2.6 
-197 
0.52,0.10 
-557, 20 




GMSB 1 


GMSB 2 


GMSB 5 


Mess. M M ess (in GeV) 
scale A (in GeV) 


10 13 
1.7- 10 5 


10 13 
1.7- 10 5 


5 • 10 14 
7.5 • 10 4 


SUSY tan % z ) 
scale (in GeV) 

A, ft 

A A ,A K (in GeV) 


8.5 
1404 
0.020,0.004 

-52,-160 


1.63 

2351 
0.50,0.43 
-446, -2300 


50 

1376 
0.010,-0.0007 
118,4645 



Table 1: Most important parameters in the mSUGRA [35[ and GMSB 36 1 benchmark scenarios 



2 ™ 2 and in case of the NMSSM 



Note, that we calculate the soft breaking masses m H , m 
in addition m? s from the tadpole equations, implying that they are not input parameters in 
CNNDecays . We do not give here the detailed soft SUSY breaking parameters but refer for 
this to [35| and 36] and/or the output of CNNDecays discussed in the appendix. 

We show the corresponding masses m, one-loop masses m 1L and particle characters C 
in Table [2J In Tables [3] and H] we present our results for the decay widths in the mSUGRA 
and GMSB scenarios, respectively. Tables [3] and H] contain beside the tree-level and one- loop 
corrected widths the correction factor 



r 1 -r c 
r° 



(103) 



which is also split in the parts 8\ = due to squark and quark corrections and 5 2 

containing the other contributions, which includes the hard photon emission for comparison 



with 12j. As already stated, the mass corrections m — > m lL presented in Table [2] are small, 
in most cases in the per-mil range. Only the light Singlino in the mSUGRA 4 scenario gets 
a large correction of 2.6% from squark and quark contributions. 

As expected the widths are larger in case the neutralino involved has either large wino 
and/or higgsino components. Note, that from ffTTl) follows that the W boson couples either 



to a wino-wino or a higgsino-higgsino combination. This explains several of at first glance 
surprising features, e.g. the fact that in mSUGRA scenarios 1 and 3 the width r(x° — > 
xtW~) is larger than r(x° — >■ W~xt) despite the smaller phase space. Another surprise 
might be the difference in 62 for the decay X2 ~* x\W~ in scenarios mSUGRA 1 and 
mSUGRA 3. However this can be understood from the differences in the Higgs sector. We 
see that the corrections are in general in the order of 1-3% but can easily go up to 10%. 
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mSUGRA 1 


mSUGRA 3 


mSUGRA 4 


GMSB 1 


GMSB 2 


GMSB 5 


Xi : m 

c 


210.79 
210.61 
B 


210.97 
21 77 

B 


89.08 
91 45 

s 


472.48 
479 38 

B 


472.53 
479 qq 

B 


203.30 
903 30 

s 


xl ■ m 

lit 

c 


387.18 

OO 1 .1U 

w 


387.47 
387 37 

OO l .o / 

w 


215.38 
91 <i ^8 

H 


620.06 
R9D Dfi 

s 


855.54 

OOO .OO 

w 


496.87 
4Qfi 81 

B 


X° ■ m 

m 1^ 
ill 

C 


971.11 
Q71 7^ 

H 


942.27 
Q41 0^ 

H 


217.09 
91 7 ^1 

H 


854.13 

8^4 4"? 

w 


2352.36 
93^9 4Q 

&00 .^O 

H 


899.60 

8QQ Q8 


xi '■ m 

c 


976.52 
975.14 
H 


943.16 
942.79 
H 


330.51 
331.01 

B 


1405.44 
1405.15 
H 


2355.92 
2354 84 
H 


1377.67 
1377.45 
H 


xi ■ m 
C 


2101.57 
2101.57 

S 


1421.70 
1421.67 
S 


608.43 
607.64 
W 


1412.41 
1411.46 
H 


4062.82 
4062.84 
S 


1383.97 
1383.04 
H 


Xi :m 
C 


387.16 
387.23 
W ± 


387.48 
387.53 
W ± 


201.36 
201.73 

H± 


854.11 
854.57 

w± 


855.53 
855.69 
W ± 


899.59 
900.14 
W ± 


xt - m 

C 


977.07 
976.69 
iT± 


947.45 
947.07 
#± 


608.40 
607.80 

w± 


1412.29 
1411.62 
#± 


2355.33 
2355.06 
iT± 


1384.24 
1383.51 
#± 



Table 2: Neutralino and chargino masses m on tree-level and m 1L on one-loop-level in GeV and main particle 
character C for the mSUGRA [35| and GMSB [36[ benchmark scenarios. 

Note, that depending on the parameters the corrections can have both signs. 

5.2. NLO corrections in case of bino and singlino like neutralinos 

It follows from (j!7p that the partial width into a jy-boson vanishes in the limit that the 
neutralino involved is either a pure bino or a pure singlino. This is the main reason that in 
Tables [3] and H] the processes containing states which are to a large extent bino or singlino 
have small widths. However, even in the limit of pure states the corresponding couplings 
are induced at the one-loop level. In this section we investigate this in more detail. We 
consider here a wide mass range and are aware that neutralinos with masses above 1 TeV 
will hardly be produced at LHC and will therefore most likely require a multi-TeV lepton 
collider such as CLIC. All the plots and decay widths are based on tree-level masses m~±o 
for neutralinos and charginos. The corresponding plots with the one-loop corrected masses 
differ only slightly and the differences are hardly visible. 
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Sc. 


Decay 


T° (in GeV) 


T 1 (in GeV) 




s 2 


#1+2 


1 


Xl -+ XiW+ 


2.153 


2.234 


1.8% 


2.0% 


3.8% 


xl^XiW + 


2.181 


2.256 


1.8% 


1.6% 


3.4% 


Xl -> XiW + 


3.206 ■ 10~ 3 


2.897- 10~ 3 


-2.6% 


-7.0% 


-9.6% 


xl -+ x*w + 


1.542 • lO" 1 


1.521 ■ lO" 1 


-1.9% 


0.5% 


-1.4% 


Xi ->■ x?^~ 


2.575 • 10" 3 


2.561 ■ 10" 3 


0.8% 


-1.3% 


-0.5% 


xi ->■ X?^" 


5.860- 10- 1 


5.766 ■ 10- 1 


0.2% 


-1.8% 


-1.6% 




2.201 


2.222 


-0.3% 


1.2% 


0.9% 


3 


xl xr^ + 


2.085 


2.153 


1.8% 


1.5% 


3.3% 




2.121 


2.181 


1.8% 


1.0% 


2.8% 




2.937- 10~ 2 


2.755 ■ 10~ 2 


-1.7% 


-4.5% 


-6.2% 


xl -+ 


1.302 


1.352 


-0.6% 


4.5% 


3.9% 


Xi ->■ x?^" 


2.951 ■ 10" 3 


2.910 ■ 10" 1 


0.7% 


-2.1% 


-1.4% 


x 2 ->■ x?^~ 


5.684- 10- 1 


5.552 ■ 10- 1 


0.2% 


-2.5% 


-2.3% 


X 2 ->■ X2^" 


2.115 


2.141 


0.6% 


0.6% 


1.2% 


4 


x2 XiW+ 


4.719 • 10~ 2 


5.080 • 10~ 2 


-0.3% 


7.9% 


7.6% 


xl -> xr^ + 


7.442 ■ lO" 1 


7.288 ■ 10" 1 


0.2% 


-2.3% 


-2.1% 


xr ->■ x?^ _ 


1.623 ■ lO" 1 


1.650 ■ 10" 1 


-0.9% 


2.5% 


1.6% 


X2" ->• X?^~ 


2.357- 10- 1 


2.291 ■ 10" 1 


0.2% 


-3.0% 


-2.8% 


X 2 ->■ X-jW" 


5.758- 10- 1 


5.586 ■ 10- 1 


0.2% 


-3.2% 


-3.0% 


X2 -»■ x° 3 w- 


6.024- lO" 1 


5.875 ■ 10- 1 


0.2% 


-2.7% 


-2.5% 


X 2 -> 


7.007- 10- 2 


6.963 • 10~ 2 


-0.3% 


-0.3% 


-0.6% 



Table 3: NLO corrections for the mSUGRA benchmark scenarios; 6 is defined in (I103p . 
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Sc. 


Decay 


T° (in GeV) 


T 1 (in GeV) 


<5i («,<?) 


s 2 


^1+2 




xl -> xr^ + 


2.891 


3.068 


-0.2% 


6.3% 


6.1% 




xl -> XiW + 


2.943 


3.114 


-0.1% 


5.9% 


5.8% 


1 


Xi -»■ x?^~ 


1.027- 10~ 2 


1.028 ■ 10- 1 


0.0% 


0.1% 


0.1% 




xr ->■ 


8.907- 10" 5 


8.942 ■ 10~ 5 


-0.1% 


0.5% 


0.4% 




X2 ->• X?^~ 


8.941 ■ 10- 1 


8.795 ■ 10- 1 


0.2% 


-1.8% 


-1.6% 




X2" ->■ X2^~ 


3.787- 10- 3 


3.697- 10" 3 


-0.5% 


-1.9% 


-2.4% 




X2~ -> x° 3 w- 


2.962 


3.126 


-0.1% 


5.6% 


5.5% 




xl -> XiW+ 


7.431 


7.241 


0.2% 


-2.8% 


-2.6% 




xl -> xr^ + 


7.442 


7.254 


0.4% 


-2.9% 


-2.5% 


2 


$ -> xr^ + 


1.164 ■ 10~ 2 


9.868 ■ 10~ 3 


-3.4% 


-11.8% 


-15.2% 




xl X2 


1.890 


1.851 


-2.0% 


-0.1% 


-2.1% 




xr ->■ xiW- 


6.150 ■ 10~ 3 


6.135 • 10~ 3 


0.0% 


-0.2% 


-0.2% 




X2 ->• x?w^- 


1.807 


1.708 


-0.7% 


-4.8% 


-5.5% 




x 2 ->■ xlw~ 


7.491 


7.235 


0.3% 


-3.7% 


-3.4% 




Xl -> XiW + 


2.283 


2.439 


-0.3% 


7.1% 


6.8% 




xl xr^ + 


2.333 


2.485 


-0.1% 


6.6% 


6.5% 


5 


xr ->■ x?^ - 


1.462 • 10~ 5 


1.453 ■ 10~ 5 


-0.3% 


-0.3% 


-0.6% 




xr -> x^" 


8.424 • 10~ 3 


8.407- 10~ 3 


-0.1% 


-0.1% 


-0.2% 




X2" ->■ XiW~ 


1.279 ■ 10~ 3 


1.242 ■ 10~ 3 


-0.6% 


-2.2% 


-2.9% 




x 2 ->■ XsiW - 


7.898 ■ 10- 1 


7.775 ■ 10- 1 


-0.1% 


-1.5% 


-1.6% 




X2~ -> xH^" 


2.339 


2.486 


-0.1% 


6.4% 


6.3% 



Table 4: NLO corrections for the GMSB benchmark scenarios; S is defined in (|103l) . 
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Figure 7: a) (left): Neutralino masses as a function of Mi and the other parameters according to mSUGRA 
1 apart from k — 0.14. The red solid line marks B whereas the other states are shown with black dashed 
lines, b) (right): Particle character for the bino state B as a function of Mi: red (solid): Bino character 
|A/ji| 2 , blue (dot-dashed): Higgsino character |A/^3 1 2 + |A/j4| 2 , black (dotted): Wino character |A^2| 2 , green 
(dashed): Singlino character |JV}s| 2 . 

5.2.1. Bino decays 

We define a state denoted by B to be bino-like if |A/}i| 2 > 0.5. We take benchmark point 
mSUGRA 1 and vary the gaugino mass M\ for our subsequent numerical investigations. 
Moreover, we shift k from 0.11 to 0.14 to disentangle different effects and to simplify the 
discussion. The corresponding neutralino mass spectrum is shown in Figure[7Ja) as a function 
of Mi and the particle character of the corresponding mainly bino-like state B is shown in 
Figure [7] b). Note that at the various crossings in plot a) the index of the corresponding 
neutralino mass eigenstate changes. 

Figure [8] a) shows the LO and NLO decay width of B — > XiW + as a function of M\. At 
Mi ~ 1 TeV the level crossing with the higgsino like states occurs giving rise to the observed 
rise of the width with Mi and the subsequent decrease. In the later case the decrease is 
strengthened by a negative interference of the higgsino and wino parts. Moreover we note 
that here and also in the corresponding figures below a Coloumb singularity occurs close to 
the kinematical threshold, e.g. close to = m-± + mw, which has to be resumed. As this 
has not been done, we start our plots slightly above this region. 

The relative size of the corrections are shown in Figure El b) where we also compare the 
third generation squark/quark contribution with the additional ones. The kinks correspond 
to the level crossings in Figure UJ Note that both parts of the correction can be of equal 
importance and that in some regions of the parameter space they cancel partly each other 
whereas there are also regions where they point in the same direction. It is remarkable that 
the loop induced corrections can be of the same order of magnitude as the tree-level widths. 
This however does not imply a break-down of perturbation theory but is a consequence that 
in the limit of a pure bino the tree-level coupling vanishes but the one-loop induced one is 
non-zero. 
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Figure 8: a) (left): LO (black, dashed) and NLO (red, solid) decay widths for B — > XiW + as a function of 
Mi for the spectrum of Figure b) (right): Correction factor S in % defined in (I103j) for B —> XiW + as 
a function of Mi: blue (dashed): Squark and quark contributions, black (dot-dashed): Other sectors, red 
(solid): Full correction. 




1200 1400 1600 1800 2000 2200 2400 1200 1400 1600 1800 2000 2200 2400 

Figure 9: a) (left): LO (black, dashed) and NLO (red, solid) decay widths for B — > x^~VF + as a function of 
Mi for the spectrum of Figure [71 b) (right): Correction factor 5 in % defined in (I103[) for B —> x^~PF + as 
a function of Mi: blue (dashed): Squark and quark contributions, black (dot-dashed): Other sectors, red 
(solid): Full correction. 

Figure [9] a) shows the widths for the decay B — > x^W^ as a function of M\. In contrast 
to the decay into Xi there is a positive interference of the wino-wino and higgsino-higgsino 
components in the LO couplings (TT7|) . Note, that the decrease of the couplings for increasing 
Mi is compensated by a factor (Mi/mw) 2 , see Equations (TT9]) -( 12^|) leading to a slight 
increase of the width with increasing M\. Also in this case the corrections can be sizable 
amounting up to about 15 percent. 
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5.2.2. Singlino decays 

The case of a singlino-like neutralino S, defined by [ iV^-g | 2 > 0.5, shows similar features 
to the case of a bino-like neutralino. However, there is one important difference: In the limit 
of a pure singlino it only couples to the doublet Higgs/higgsino states and the singlet Higgs 
boson. Hence one expects that the squark/quark contributions to be of less importance 
compared to the bino case. 




Figure 10: a) (left): Neutralino masses as a function of k and the other parameters according to mSUGRA 
4 apart from A = 0.01. The green solid line marks S whereas the other states are shown with black dashed 
lines, b) (right): Particle character for the singlino state S as a function of k: red (solid): Bino character 
|Aji| 2 , blue (dot-dashed): Higgsino character |^Vj 3 1 2 + |A/^4 1 2 , black (dotted): Wino character |A/^2| 2 , green 
(dashed): Singlino character \Afj5\ 2 . 




Figure 11: a) (left): LO (black, dashed) and NLO (red, solid) decay widths for S — > XiW + as a function 
of k for the spectrum of Figure [TU1 b) (right): Correction factor S in % defined in (I103P for S — > XiW + 
as a function of k: blue (dashed): Squark and quark contributions, black (dot-dashed): Other sectors, red 
(solid): Full correction. 
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Figure 12: a) (left) LO (black, dashed) and NLO (red, solid) decay widths for S —> x^W + as a function 
of k for the spectrum of Figure [TU1 b) (right): Correction factor 5 in % defined in (|103[) for S — > xJW^ 1- 
as a function of k: blue (dashed): Squark and quark contributions, black (dot-dashed): Other sectors, red 
(solid): Full correction. 

For the numerical investigation we choose the benchmark scenario mSUGRA 4, but 
reduce A to A = 0.01 to increase the singlino character. We vary k between 2 • 10~ 3 and 
6 • 10 -2 leading to singlino masses between 100 GeV and 2.5 TeV. In this rather light particle 
spectrum the higgsino-like chargino has a mass of m~± =201 GeV and the wino-like chargino 
has a mass of m-± = 608 GeV. The mass spectrum of the neutralinos is given in Figure fTOl a). 

In Figure [TT] we show the details for the decay S — > XiW + . Note that despite the 
decrease of the coupling due to the decrease in the wino and higgsino components with 
increasing k one gets a net increase of the width due to the (m^/m^) 2 factor. In the 
right plot one sees clearly that contributions of third generation squarks and quarks are less 
important compared to the bino case. However, the remaining contribution can amount up 
to about 10 percent. Also the decay into the heavier chargino shows similar features as can 
be seen from Figure [121 The main difference is that the threshold effects due to on-shell 
intermediate states are more pronounced in this case and are mainly caused by sleptons and 
Higgs bosons at irtg « 1 TeV and by squarks at m§ w 1.6 TeV. 

5.2.3. Chargino decays 

In this section we present the decays of the wino-like chargino W + in a bino- or singlino- 
like neutralino Xi2> being the two lightest neutralinos in the benchmark scenario GMSB 5. 
We depart from the original parameters by setting Mi = 300 GeV and \i = 600 GeV to get 
a rather light particle spectrum. We vary the gaugino mass M 2 between 100 and 2000 GeV. 
The neutralino mass spectrum is given in Figure [JJJ] a) and the chargino mass spectrum in 
Figure [TBI b). The two light neutralinos have a mass of m^o = 89 GeV for the singlino-like 

neutralino and m^o = 298 GeV for the bino-like neutralino. 

2 

In Figures [T41 and [TBI we show the details of the decays W + — >■ x?2^ + - Note that the 
peaks close to M 2 ~ 620 GeV are due to the level crossing of the wino-like states with the 
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higgsino-like states. The overall features of the widths and corrections are similar to the 
case of neutralino decays. In general the corrections are in the order of a few percent except 
for a region close to M<i = 1.15 TeV for the decay W + — > X2^ + m Figure IT5l where the 
tree-level couplings to x® near ly vanish due to a negative interference between the wino and 
higgsino contributions. We find that the third generation squark/quark contributions are in 
general somewhat smaller than the remaining ones. 




Figure 13: a) (left): Neutralino masses and b) (right): Chargino masses as a function of M%. The other 
parameters are as GMSB 5 apart from M\ — 300 GeV and /i = 600 GcV. The red lines in b) correspond to 
the wino like states and the two blue ones in a) to the singlino state Xi an d bino state xl- 
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Figure 14: a) (left): LO (black, dashed) and NLO (red, solid) decay widths for W + — > XiW + as a function 
of Mi for the spectrum of Figure [TJl b) (right): Correction factor 5 in % defined in f|103[> for W + — > x\W + 
as a function of M^: blue (dashed): Squark and quark contributions, black (dot-dashed): Other sectors, red 
(solid): Full correction. 
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Figure 15: a) (left): LO (black, dashed) and NLO (red, solid) decay widths for W + — > X2^ + as a function 
of M2 for the spectrum of Figure [T31 b) (right): Correction factor S in % defined in (I103|) for W + — > x§VF + 
as a function of M^- blue (dashed): Squark and quark contributions, black (dot-dashed): Other sectors, red 
(solid): Full correction. 

6. Conclusion 

In this paper we presented the complete electroweak one-loop corrections to the partial 
widths for two-body decays of a chargino (neutralino) into a VT-boson and a neutralino 
(chargino). The calculations were done for the MSSM and the NMSSM using an on-shell 
scheme and we have used the so-called pinch technique to achieve gauge invariance. Moreover 
we based the calculation of one-loop masses for neutralinos and charginos on a reasonable 
renormalization of the physical parameters, which are present on tree-level. 

Whereas the mass corrections to neutralinos and charginos are typically in the per-mil 
range, the corrections to the decays are in the order of 1-10 percent, in particular if the 
involved neutralino is either wino or higgsino like. However in case of a small leading order 
decay width the corrections can be 50 percent or larger. This is the case if either there is 
a negative interference between the wino and the higgsino part of the couplings or if the 
neutralino is either mainly bino- or singlino-like, because for a pure bino or a pure singlino 
the tree-level width is identical zero and can only be induced at the 1-loop level. Last but 
not least we provide the publicly available program CNNDecays for the numerical evaluation. 
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Appendix A. Formulas: Vertex corrections 



In this section we present the generic formulas for the next-to leading order vertex con- 
tributions for the 't Hooft-Feynman gauge £y = 1 (V = W, Z) . The formulas for ^ 1 are 



contained in the program CNNDecays 34]. The formulas for the self-energies can be found 



in 



38[, from which also the derivatives with respect to p 2 can be calculated. They are also 
included in CNNDecays [j3 |. 

Below we present the formulas for the generic contributions to the matrix element My 
shown in Figure HI In addition we give the particle combinations to be inserted in these 
diagrams for the decay x° X~W + neglecting generation indices. We use the following 
notation: H stands for one of the 2(3) scalar Higgs bosons, A for one of the 2(3) pseudoscalar 
Higgs bosons including the Goldstone boson and H ± for one of the 2 charged scalars including 
the Goldstone boson. The indices of couplings and masses in the generic formulas have to 
be understood in the following form: and F denote the decaying and outgoing fermion, 
W the external PF-boson, whereas F, S, Si^, V or V X;2 represent possible internal 

fermionic, scalar or vector particles. It is understood implicitly that one has to sum over 
possible flavour and generation indices of the internal particles. All contributions have the 
same generic structure: 

My = T ^u( Pl )r(PLM 1 + P R M 2 )u(k)e;(p 2 ) (A.l) 
u( Pl ) (P L Mg + P R M£) u(k)e*(p 2 ) 
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We start with the vertex in Figure H] a) containing two internal fermions and one scalar. 
We use the following abbreviations for the couplings: 

Oi = O ffv>l (F 2 , F u W), 2 = F fv,r(F 2 , Ft, W) (A.2) 
3 = FFS , L (Fu Fi, S), 4 = FFS , R (Fu F u S) (A3) 
5 = FFS , L (Fo, F 2 , S), O e = O ffs ,r(Foi F 2) S) (A.4) 

In the formulas below the Passarino-Veltman integrals have as arguments: Bq{w^ : m F2 ,m Fl ) 
and Cj, Cjj(TO^, m 2 , m 2 , m F2 , m F , rn 2 s ). Possible particle insertions in the notation SFiF 2 
are given by Hx ^^ Ax^X^i H ± x ± x°, vvl, llu, uud, ddu. We get: 

M lta =-0 1 [O 3 m Fl (O 6 m F2 C - O h m C x ) + O 4 O G m F2 mi(C + C x )] (A.5) 
+ 1 5 m C 2 {0 3 m Fl + O^rrii) + 2 [0 3 {0 6 {m 2 C Q + m 2 C x + B - 2C 00 
-m 2 o (C + C X + C 2 )) - O 5 m F2 m o (C + C x + C 2 )) + O^O^m^rriid] 

M %a =O x {OsO.mp^C, + 4 [O 5 (m 2 C + m\C x - m 2 o (C + C X + C 2 ) (A.6) 
+B - 2C 00 ) - O 6 m F2 m o (C + C X + C 2 )]} 
+ 2 {0 3 mi{0 6 m C 2 - O 5 m F2 (C + C x )) 
+0 A m Fl (0 6 m (Ci + C 2 ) - O 5 m F . 2 C )} 

M^ a =pt2 {O x 4 [0 5 m (C 12 + C 2 + C 22 ) - O 6 m F2 (C + C X + C 2 )} (A.7) 

29 



+0 2 6 (0 4 m Fl Ci - 3 miC 12 )} 

+ p%2 {0 2 % [0 3 m,(d + C 12 ) + Ownd] - 1 O^O i m {C 1 + 2C 12 )} 
M£ a =p^2 {0 1 5 (0 3 m i r 1 C 1 - 4 mid 2 ) (A.8) 
+0 2 3 [0 6 m (d 2 + d + C 22 ) - O 5 m F2 (C + d + C 2 )]} 
+ p£2 {d0 5 (0 3 m Fl d + 4 m;(d + C 12 )) - 2 O^O^m {C x + 2C 12 )} 

For the vertex in Figure H] b) we define: 

1 = O ssv (S 1 ,S 2 ,V) (A.9) 

2 = O ffs ,l{F, F h Si), 3 = O FFSiR {F, F u SO (A.10) 
d = FFS , L (F , F, S 2 ), 5 = O ffs ,r{F 0) F, S 2 ) (A.ll) 

For completeness we note that in case of 0\ the following internal momentum combina- 
tion appears {p$ 1 — ps 2 ) over which of course has been integrated. The Passarino-Veltman 
integrals below have the arguments Cj, Cij{m 2 , m^, rn 2 , m 2 F , rrvg 2 ,rrig ). Possible particle in- 
sertions in the notation FSiS 2 are given by x Q HH ± , x ± H ± H, x°AH , x H A, vvl, UP, 
uud, ddu. We obtain for the different parts: 

M 1>b =2O 1 O 2 O 5 C 0Q (A.12) 

M 2tb =20,030,000 (A.13) 

K,b =^d {O 3 (O 5 m F (C + 2C 2 ) - 4 m c (d + 2C l2 )) (A.14) 

-0 2 5mi (C 2 + 2C 22 )} - 2p»0 1 {0 2 5 m t (C l2 + C 2 + C 22 ) 

+0 3 [0 4 m (d + di + C 12 ) - O 5 m F (C + d + C 2 )]} 

M£ b =2pt0 1 {0 2 [O 4 m F (C + d + C 2 ) - 5 m (d + C n + C 12 )} (A.15) 

-0 3 4mi (C 12 + d + C 22 )} + P^O, {0 2 [O 4 m F (C + 2C 2 ) 

-0 5 m o (d + 2C 12 )] - 3 A m t {C 2 + 2C 22 )} 

For the vertex in Figure H] c) we define: 

Ox = O ffv>l (Fi, F { , V), 2 = O ffv>r (Fi, F h V) (A.16) 
3 = O mL (F 2 , F u W), 4 = O FFViR (F 2 , Fi, W) (A.17) 
5 = FFV , L (F , F 2 , V), 6 = O ffv>r (F , F 2 , V) (A.18) 

The Passarino-Veltman integrals below have as arguments d, ^(mfy, m 2 , m 2 , m F ,m F ,my) 
and Bo{m 2 v , mi, m F ). Possible particle insertions in the notation VFiF 2 are given by 
Zx ^, W^x . We get: 

M ljC =20i0 5 {0 3 [m^C + m?(d - d) - ro*(C + d + 2d) (A.19) 
+m^/d + -Bo - 2do] - 4 m Fl m F2 Co] - 20 2 A 6 mim C 2 

M 2iC =20 2 6 {0 4 [m^C + m, 2 (d - d) - m*(Co + d + 2d) (A.20) 
+m^d + Bo - 2do] - 3 m fl m Fi C } - 20 1 3 5 m i m C 2 
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=4^ {0 2 [0 3 5 m Fl C 2 + 4 (0 5 m F2 C 2 + 6 m (C 2 + C 22 ))} (A.21) 
-dOsCWCs + C 12 )} - 4p% {0 2 [0 3 5 m Fl d + 4 (0 5 m F2 (Co + C x ) 
+0 6 m (C 1 + C n + C 12 + C 2 ))] - O^O.m^d + C n )} 

M£ c =4p? {Ox [0 3 (Or m (C 12 + C 22 ) + 6 m F2 C 2 ) + 4 6 m Fl C 2 } (A.22) 
-0 2 0,0 6mi (C 2 + C 12 )} - 4p£ {Oi [0 3 (0 5 m (C 1 + C n + C 12 + C 2 ) 
+O 6 m F2 (C + d)) + O 4 6 m Fl Ci] - O^C^d + C u )} 

For the vertex in Figure H] d) we define: 

Oi = O FF v,l{Foi F, V), 2 = FFV , R (F 0} F, V) (A.23) 

3 = Offs,l{F, F h S), A = O ffs , r (F, F h S) (A.24) 

5 = O svv (S,W,V) (A.25) 

The Passarino-Veltman integrals below have as arguments Cj, Cyim^,, rriyy, ™ F i 
Possible particle insertions in the notation .FSV are given by x ± H ±/ y, x ± H ± Z, x°HW, 
X°AW. We get: 

M M =O 5 (O 1 O 3 m F C - O^m.C, + O-An^d) (A.26) 

M 2id =O 5 (0iO \m C x - O^m.C, + O 2 O 4 m F C ) (A.27) 

M 3 M d =2rfO 1 O 4 5 C 1 (A.28) 

M£ d =2^0 2 3 5 C a (A.29) 

For the vertex in Figure H]e) we define: 

Oi = O FF v,l{F, Fi, V), 2 = O ffv , r (F, F h V) (A.30) 

3 = O ffs ,l(Fo, F, S), 4 = FFS , R (F , F, S) (A.31) 

5 = S vv(S,W,V) (A.32) 

The Passarino-Veltman integrals, which appear in the following formulas, have as arguments 
Ci, Cij(mf, m^y, ml, m F , rriy, m 2 s ). Possible particle insertions in the notation FVS are given 
by jfZH*, x^WH, x±WA. We obtain: 

M li6 =0^-0x0^^ + O 1 O 4 m F C + O^^d) (A.33) 

M 2 , e =0^0x0^^ + O 2 O 3 m F C - 2 A m C 2 ) (A.34) 

Ml e =2(rf + p%)0 2 A 5 d (A.35) 

M£ e =2(p? + p£)Oi0 3 B d (A.36) 

For the vertex in Figure H]f) we define: 

1 = O FF v^l{F, Fi, V,), 2 = O FF v.r{F, F h V x ) (A.37) 

3 = O ffv ,l(Fo, F, V 2 ), 4 = O ffv , r (F , F, V 2 ) (A.38) 

o 5 = o wy (w,y 1 ,y 2 ) (A.39) 

31 



For completeness we note, that O5 has the following internal momentum contribution over 
which has been integrated: ((py. — Pw)g ucr + •••)• The Passarino-Veltman integrals, which 
appear in the following formulas, have the arguments Cj, Cij(m 2 0) m^, wi? , Trip, Triy ,my ) and 
Bo(myy, m y 2 i m y i)- Possible particle insertions in the notation FV1V2 are given by X W^f, 
X°ZW, x±WZ. We obtain: 

Afx,/ =0 5 {Ox [O 3 (2m|C + m?(2(7 1 + 3C 2 ) (A.40) 

+m2(3d + 2C 2 ) - 2m 2 w {C 1 + C 2 ) + 2B + 4C 00 ) + 3O 4 m F m o C ] 

+30 2 mi [O 3 m F C + 4 m (d + C 2 )]} 
M 2i/ =0 5 {30im i (0 3 m (C 1 + C 2 ) + O 4 m F C ) + 2 [3CWm o C (A.41) 

+ O 4 (2m 2 F C + m?(2d + 3C 2 ) 

+m*(3C 1 + 2C 2 ) - 2m 2 w/ (C 1 + C 2 ) + 25 + 4C 00 )] } 
M£ f = - 20 5 K [0 1 3 m l (2(C 12 + C 22 ) - C x ) + 2 (30 3 m F (d + C 2 ) (A.42) 

+0 4 m (2(C 11 + C12) - C 2 ))\ + p% [Ox0 3mi (C 2 + 2C 22 ) 

-0 2 (0 4 m (Cx - 2C 12 + C 2 ) - 30 3 m F C 2 )]} 
M£ f = - 20, K [Ox(0 3 m (2(Cxx + C 12 ) - C 2 ) + 3CW(Ci + C 2 )) (A.43) 

+0 2 4 m i (2(C 12 + C 22 ) - Cx)] + p£ [0 2 4 rrii(C 2 + 2C 22 ) 

-Ox(0 3 m {Cx - 2C 12 + C 2 ) - 30 4 m F C 2 )]} 



Appendix B. Hard photon emission 

Here we present the formulas for the hard photon emission for the Feynman graphs given 
in Figures [5] and El We use the notation of [14| for the Bremsstrahlung integrals, which are 
defined for the decay of a particle with mass mo and momentum po into two particles with 
masses mi and m 2 and momenta p\ and p 2 and a photon with momentum q by 

rij, \ 1 f d 3 p l d 3 p 2 d 3 q (4) (±2p i q)(±2p j q) 

where the minus signs refer to the momentum po of the initial particle. In this context we 
need the integrals J^(mj, m , m w ), allowing us to write the final result in the form 



327r 2 m i m 2 v 



Owl? + \0 WR \ 2 ) fix + (0 WL 0* WR + 0* WL WR ) Q 2 ] , (B.2) 



where we have introduced abbreviations 



Qi = Q 2 flxu + Ql^iioo (B-3) 

^2 — Ql^2ii + Ql^2oo 
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with the charges Qi and Q Q of in- and outgoing particles being either or 1 depending on 
the process considered. Due to the presence of left- and right-handed couplings to the W 
boson in ( TT6|) the final result of the three-body does not factorize in the two-body width 
times a corrections factor. The different parts are given by: 



Vt m =21 (m 2 + m 2 + 2m 2 w ) - 4 [m 2 w (m 2 + w 2 ) + (m 2 - m 2 f - 2m A w ] 

■ [l + I 00 m 2 + m 2 w (I 02 + I22) + io2(m* - m )(nii + m ) + I 2 ] 
+ 2I 2 (m 2 + m 2 + 2m 2 w ) - 8m 2 w (l!g + l\) 
n loo = - 8Im 2 w - Am 2 w [m 2 o (m 2 (2I 01 + 2I 02 



n 



2ii 



n 



2oo 



+ m (-2/ 01 - 2/ 02 + I 



'22; 



h In - 2/ 22 ) + A + J 2 ) 
m?(/i + I 2 + /22m, 2 ) + 7° + 21 



(B.4) 



(B.5) 



- J 22 ) + 4m^ [-m 2 (7 i + J 02 
/02) - / 22 ) + 2/i + 2/ 2 ] 



+4/ 2 ° + 2/2 1 ] 

- 8m^(/ i + /02 

+m 2 (3(/ i + 

- 4/oim^ + 12/oim- m 2 , - 12/ im 2 m^ + 4J im^ 
+ 12/ 02 mfm 2 - 12/ 02 m 2 m^ + AI 02 m 6 o - 4J 1 m l 4 

A A O 0/1 

+ %l\m i m — AIim — AI\\m i m + 8Jnm i m 

- 4J u m^ - 2J 1 °m 2 - 2/ 1 °m 2 - 4I 2 m* + M 2 m 2 m, 



2/n + / 22 ) 
4/02^6 



4J 2 m„ 



--2m imo [-21 + I2m 2 w (l + m 2 (I 00 + I 02 ) - I 02 m 2 o + I 2 ) 

+I2m 4 w (l 02 + J 22 ) - 2J 2 ] 
-Ammo [6m^(-mo(/ i + -^02 - Ai) + m 2 (/ i + -M + A + ^2) 

-6m^(/ i + /02 - J22) + i?] 



This result was obtained by using FeynArts [39[ and FormCalc 4^ . 



01 

22 



(B.6) 
(B.7) 



Appendix C. One- and Two-point functions 



Next we provide some useful formulas for the one- and two-point loop functions 41 



As only part of them can be found in different places in the literature we collect here the 
complete set, in particular the derivatives of the two-point functions £?ooi and 5m. We 
remind that the UV divergent part of the integrals can be expressed as 



1 



A = 7 + ln(47r) 



(CI) 



where 7 ~ 0.577 denotes Euler's constant. For the notation of the Passarino-Veltman 



integrals we follow [40|, |42 



Appendix C.l. Scalar integrals 

The one-point function Aq is given by 



A n (m 2 ) 



m 



A 



I ^\— 2 



(C.2) 
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and the two-point function B by 



(p 2 , ml ml) = A + 2 + In ( + ^1 In ( ^) - ^ ( 1 - r) In r, (C.3) 

\m\m2 J p \m± j p z \r J 



where r and - denote the negative roots of the polynomial 



2 1 2 2 
2 mf + mi — p L 

x H x + 1 — (x + r)\x + 

m\m2 



The derivative of B with respect to p 2 yields 

dp 2 



B (p 2 ,m 2 1 ,m 2 2 ) := -^B (p 2 ,ml,m 2 2 ) 



2 2 
mj — wij 



ln[^) + 

mi 



m\m2 ( 1 



r ) lnr - ( 1 + 



r 2 + x 



lnr 



(C.4) 



Appendix C.2. Tensor integrals 

Lorentz covariance in d dimensions allows to decompose the tensor integrals in terms of 
scalar integrals which are given by: 



I 1 

A 00 (m 2 ) = -m 2 A (m 2 ) + -m 4 
4 o 

1 



B 1 (p 2 ,m 2 1 ,m 2 2 ) = — [( 



m — m 



j)(B (p 2 ,mlm 2 2 ) - B (0,mlm 2 2 ))] 



-B (p 2 ,ml,m 2 2 ) 



B 00 (p 2 ,ml,m 2 2 ) 



A (ml) + (p 2 -m\ + mj)B 1 (p 2 ,m 2 ,m 2 .) 



+2mlB (p 2 , ml, ml) + ml + ml - -p 2 



1 

3p 2 



B n (p ,m 1: m 2 ) = — A (m 2 ) - m l B Q (p ,m 1 ,m 2 ) 



-2(p 2 — m 2 , + ml)Bi(p 2 , m\, m 2 2 ) + -(p 2 — 3ml — 3m 



Booi(p 2 ,ml,ml) = - 2m 2 1 B 1 (p 2 ,m\,m\) - A (m 2 2 ) 



+ {p 2 -m 2 2 + ml)B n {p 2 , ml, m 2 2 ) - -(2m? + 4m^ - p 2 ) 



B in (p 2 ,ml,ml) 



1 



Ap 2 



A (ml) + 3(p 2 - ml + ml)B u (p 2 , ml, ml) 



(C.5) 
(C.6) 

(C.7) 
(C.8) 
(C.9) 
(CIO) 



+2mlB 1 (p 2 ,m 2 1 ,ml) - -(2ml + 4m 2 _ P 2 ) 
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PV integral 


UV behaviour 


PV integral 


UV behaviour 


A 


m 2 A 


^00 


±m 4 A 


Bo 


A 


B 1 


-1A 


Boo 


^(3m 2 + 3m 2 - p 2 )A 


B n 


|A 


Boot 


±(-2m 2 -Am 2 + P 2 )& 


Bm 


-|A 


Coo 


|A 


Cooi 




Cq02 








-Boo 




Bqoi 


^A 



Tabic C.5: UV divergent parts of the Passarino-Veltman integrals. 



Appendix C.3. Special cases for B functions 

The following special cases turn out to be useful in the numerical evaluation. Here we 
give only the finite parts and summarize the UV divergent parts of the functions appearing 
in the calculation in Table IC.5I 



£>o(0, m 2 , m 2 



£ (p 2 ,0,m 2 ) 
B (p 2 , m 2 , m 2 ) 

D I 2 2 2^ 

Bo{m ,m ,m ) 



Bo(0,0,m 2 )=B (Q,m 2 ,0) = l+\n x „ 

m z 

. 1 



2 2 

m\ — m 2 



m i 111 ( 32 l _ m 2 111 f 32 



2 



,00(0, m 2 ,m 2 ) = In 



Q 2 



S o (^,0,0) = 2 + ln ^ + m 



Q 2 



v 



B o (p 2 ,m 2 ,0) = 2 + In 



Q 2 



+ 



2 2 
m — p 



In 1 



P 



2 -flu ( ^ -^(--r )\nr 
171 J p \r 

2 + In. I ^4 I - " 



(C.ll) 
(C.12) 
(C.13) 
(C.14) 
(C.15) 
(C.16) 
(C.17) 



Appendix C.4- Derivatives of the B functions 

First we present the general results for the derivatives and then the special cases. As 
above we give here only the finite parts as the UV divergent parts are given in Table IC.5I 
Bo(p 2 ,mf,ml) is given in (IC4[) . 

1 



B x {p ,m 1 ,m 2 ) 



2p 4 



[(ml - m 2 2 )Bo(p 2 , mj, ml) + (ml - ml)B (0, m 2 , ml) 
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2/2 2 , 2\ r> I 2 2 2\ 

— p [m 1 — m 2 + p )Bo(p , m 1 ,m 2 ) 
1 



(C.18) 



B 00 (p 2 ,ml,ml) [-3(m 2 - m^^o^m^m 2 ,) 

+ 3(m 4 - 2m?m^ + - p 4 )S (p 2 , mj, m\ 
-P 2 (3k(p 2 , m\, m 2 2 )B (p 2 , m 2 , ml) + 2p 2 ) 

B u (p 2 ,ml,ml) [2(mi - m 2 )(m 1 + m 2 )(2m\ - 2m 2 +p 2 )S (0, m 2 , m 2 ) 



(C.19) 



- 2(p 2 (m 2 - 2m 2 ) + 2(m 2 - m 2 2 ) 2 )B (p 2 , m 2 , m 2 ; 

+ 2p 2 {p 2 {m\ - 2ml) + ( m i - m 2) 2 + P 4 ) B o(p 2 , m\, ml) 
-2p 2 A Q (ml) + p 2 {m\ + m\ 



(C.20) 



S 00 i(p 2 , m\, m 2 2 ) 



144 p e 



[6(m 2 -m 2 )(2(m 2 -m 2 ) 2 -p 2 (m 2 + 2m 2 )) J B (0,m 2 ,m 2 ; 



B in (p 2 , m\,m 



+ 6(p 2 (m 2 - mjiK + 3m 2 .) - 2(m 2 - m 2 ) 3 + p 6 )B {p 2 , m\,m 2 2 ) 
+ 6p 2 (m 2 1 -m 2 2 +p 2 ) 

{-2p 2 {m\ + m 2 ) + (m 2 - m 2 ) 2 + p 4 )B (p 2 , mj, m 2 ) 
+6p 2 (m 2 - m 2 )A (m 2 ) + p 2 (3m 4 - 3m 4 + 4p 4 )] (C.21) 

[3 (m 2 - m 2 2 ) (2p 2 (mj - 2m 2 . ) 



I2p 



+ 3(m 2 - m 2 2 ) 2 + p 4 )5 (0, m 2 , mi 



— 3(p 4 (m 2 — 3m 2 ,) + 2p 2 (m 2 — Sm^mJ — ml) 

+ 3(m 2 l -m 2 2 f)B {p 2 ,m 2 1 ,m 2 2 ) 
+ 3p 2 (m 2 - m 2 +p 2 )((m 2 - m 2 ) 2 - 2m 2 p 2 + p 4 )B (p 2 , m\, m 2 ) 
-6p 2 A (m 2 )(m 2 -m 2 +p 2 ) 
+p 2 (3m 4 + 2m 2 p 2 - 3m2 + km 2 2 p 2 )\ 

In subsequent formulas for p 2 = we use the abbreviations 



(C.22) 



K x = log 
and get for m 1 ^ m 2 ^ 
S (0,m 2 ,m 2 ) 
5i(0,m 2 ,m|) 
5oo(0,m 2 ,m^) 



Q 1 



mi 



'^2 = log ( ^ ) , K 3 = log 



1 



If 



2{m\ — m 
1 

6(mf — m 2 ,) 4 
1 



m^ 



, 2m 2 m 2 Ks 



mi 



ml 



m x + m 2 + 



mi 



mi 



\-3m\ml [2K 3 + 1) - 2m\ + 6m 2 m 4 - m|] 



[- (6Kx + 5) m\ + 9 (27^ + 3) 



72(mf -m|) 3 
-9 (2K 2 + 3) m 2 m 4 + (QK 2 + 5) m°] 
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m x mr z 



(C.23) 

(C.24) 
(C.25) 

(C.26) 



J3n(0, m?,m|) 



[Am\ml (6K 3 + 5) 



Sooi(0,m?,m|) 



Siii(0,m?,m|) 



24(m? - m 2 ,) 5 
+6m^ - 36m^m2 + \2m\m\ - 2ml] 

288W-^)« [(12A '' + 13 )-5-8(6if, + 11) mV? 
+ 36 (2^2 + 3) m£m£ - 8 (6# 2 + 5) m 2 ™^ 
+ (12K 2 + 7) mjj] 

1 [-hm\m 2 2 {12K 3 + 13) 



60(m 2 - m 2 ,) 6 
-12m}° - \2Qm\m\ - mm\m G 2 + 20m?ml - 3m^°] 

and for the remaining cases: 

1 



B (0, m 2 ,m 2 ) 
B (0, 0,m 2 ) = 



6m 2 



S (0,w 2 ,0) 



2m 2 
1 



2m 2 



B 00 (0,m ,m 
B 00 (0,0,m 2 ) = — 
B 00 {0,m 2 ,0) 



12 



m z 



1 

72 
1 

72 



6 log 




1+5 








6 log 













S 00 i(0,0,m 2 ) 



24 



288 



12 log 



Q 2 



rrr 



+ 7 



5ooi(0,m 2 ,0) = — 



121og(^)+13 



Bi(0, 0,m 2 ) = 
l?i(0,m 2 ,0) = 
fin(0, m 2 ,m 2 ) 

^(0,0, m 2 ) = 
B u (0,m 2 ,0) = 
Sm(0, m 2 ,m 2 
flm(0,0,m 2 ) = 
fi m (0,m 2 ,0) = 



1 



12m 2 



6m 2 



3m 2 



1 



20m 2 
1 



12m 2 
1 



Am 2 



1 



30m 2 
1 



20m 2 
1 

5m 2 
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C.28) 
C.29) 

C.30) 
C.31) 
C.32) 
C.33) 

C.34) 

C.35) 

C.36) 

C.37) 
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Appendix D. The program CNNDecays 



We provide for the numerical evaluation of the NLO corrections to the decays xt ~~ 



XiW and Xi ~~ >* XkW the program CNNDecays which can be obtained from 34j. It is 
written in Fortran 95 and based on SPheno 43J. The program folder contains the following 
sub-folders: 

• callcorrections: routines to combine the generic routines contained in corrections 
with the model dependent information concerning masses and couplings. 
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corrections: generic NLO routines, which are provided in i?£-gauge and 't Hooft- 
Feynman gauge, as well as the loop functions which are not contained in the SPheno 
package. 

couplings: NMSSM couplings which have been generated by a Mathematica code 



and then were cross-checked with the program Sarah [44 



• oneloop: contains the main program CNNDecays . f 90 and the main module for the cal- 
culation Renormbasic . f 90. The later one can also be used to implement the package 
in other programs. The calculation of wave-function renormalization constants and 
counterterms is included in wavemassrenorm. f 90. The module Bremsstrahlung.f 90 
contains the calculation of the hard photon emission. 

• sphenooriginal: necessary parts of SPheno. 

Before compiling it might be necessary to adjust the f90-compiler and the corresponding 
flags in the Makefile which is placed in the main folder. Using make the program CNNDecays 
is created which is stored in the sub-folder bin. 



The input and output files are based on the SUSY Les Houches Accord (SLHA) [28|, [29 . 
Concerning the input, which is expected to be given at the electroweak scale, there are two 
main differences with respect to the SLHA: 

1. The entries of the block EXTPAR are interpreted as effective on-shell values for the 
masses and mixing entries. Therefore the entry setting the scale is ignored. 

2. A new block called NLOPAR has been created containing the information to check for 
the gauge and renormalization scale independence of the results. The program allows 
to use only the UV divergent parts of the Passarino-Veltman integrals. Moreover the 
divergence itself can be set to an arbitrary value. By varying the photon mass, one 
can simply check the IR finiteness. In addition the gauge parameter £y can be set 
to an arbitrary value and it can be chosen, whether i?g-gauge or 't Hooft-Feynman 
gauge should be used for the photon, the Z-boson and the W-bosoia independently. 
Note, that the renormalization scale Q in NLOPAR only affects the scale within the 
Passarino-Veltman integrals and does not imply any running of the parameters of the 
block EXTPAR. Last but not least, one can choose whether LO or NLO neutralino and 
chargino on-shell masses are used for the calculation of the processes, meaning they 
enter as external as well as internal masses. 

We use for both models, the MSSM and the NMSSM, the entry 23 in the block EXTPAR to 
provide /x. In case of the NMSSM it is used together with entry 61 to calculate the singlet 
vacuum expectation value v$- Below we give an example input file LesHouches . in based 
on the benchmark scenario mSUGRA 1: 



Block MODSEL 

3 1 
Block SMINPUTS 

1 1.27920000E+02 



# Select model 

# NMSSM 

# Standard Model inputs 
# ALPHA.EM" -l(MZ) 
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2 1.16639000E-05 

3 1.17200000E-01 

4 9.11870000E+01 

5 4.21400000E+00 

6 1.71400000E+02 

7 1.77700000E+00 
Block MTNPAR 

3 10 

4 1 
BLOCK EXTPAR 

1 2.11635141E+02 

2 3.91898115E+02 

3 1.11230823E+03 

11 -1.42395369E+03 

12 -2.61046378E+03 

13 -1.77741018E+03 
23 9.68523016E+02 

31 3.77391179E+02 

32 3.77391179E+02 

33 3.65780798E+02 

34 2.57517852E+02 

35 2.57517852E+02 

36 2.21138594E+02 

41 1.02413355E+03 

42 1.02413355E+03 

43 8.46048325E+02 

44 9.86146013E+02 

45 9.86146013E+02 

46 5.65558510E+02 

47 9.81632542E+02 

48 9.81632542E+02 

49 9.66133394E+02 

61 1.00000000E-01 

62 1.08485437E-01 

63 -9.59966990E+02 

64 -1.58051889E+00 



BLOCK NLOPAR 



1 

2 0.00000000E+00 

3 1. 

4 9. 

5 1. 

6 1 

7 1 

8 1 

9 



00000000E-05 
11870000E+01 
00000000E+05 



# GF 

# ALPHAJ3 (MZ) 

# MZ 

# MB(MB) 

# MTOP (POLE MASS) 

# MTAU 

# Input parameters 

# tanb 

# sign (mu) 

# Ml 

# M2 

# M3 

# ATOP 

# ABOT 

# ATAU 

# MU 

# M_eL 

# M_muL 

# M.tauL 

# M.eR 

# MjnuR 

# M.tauR 

# M_qlL 

# M_q2L 

# M_q3L 

# M.uR 

# M.cR 

# M.tR 

# MjdR 

# M_sR 

# M_bR 

# LAMBDA 

# KAPPA/ 2 

# TJjAMBDA/LAMBDA = AXAMBDA 

# TJ<APPA/KAPPA = AJ<APPA 

# UV divergence : 1 = only UV div . parts 

# UV divergence : Delta 

# IR divergence : Photon regulator mass 

# Renormalization scale : Q for NLO calc . 

# Gauge dependence : Xi 

# Xi = 1 for photon , otherwise set to 

# Xi = 1 for Z, otherwise set to 

# Xi = 1 for W"\pm, otherwise set to 

# NLO masses for process : uses LO masses 



1 uses NLO masses 



A successful run creates the output file CNNDecays .dec. In the example we give only the 
crucial information. We store in the SLHA block MASS the NLO masses of neutralinos and 
charginos, whereas the LO masses are only part of the screen output. In the SLHA block 
DECAYTREE the LO decay widths T in GeV are shown. The corresponding NLO decay width 
T 1 in GeV are given in the SLHA block DECAY: 



# SUSY Les Houches Accord 2 — Neutralino + Chargino NLO Decays into W boson 

# S. Liebler , unpublished 

# in case of problems send email to sliebler@physik . uni— wuerzburg . de 

# Created: 02.03.2011, 19:14 

Block SPINFO # Program information 

1 CNNDecays # spectrum calculator 

Block SMLNPUTS # SM parameters 
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Block MODSEL # Model selection 

3 1 # NMSSM 

Block MINPAR # Input parameters 

3 1.00000000E+01 # tanb at 

Block EXTPAR # 



m_Z 



Block gauge 



Block MASS # Mass spectrum 
# PDG code mass 



1000022 
1000023 
1000025 
1000035 
1000045 



2.10611052E+02 
3.87101530E+02 
-9.71754981E+02 
9.75135414E+02 
2.10157308E+03 



particle 

# chiOl 

# chi02 

# chi03 

# chi04 

# S 



1000024 
1000037 
DECAYTREE 

# BR 
9.79638284E-01 
2.03617164E-02 

DECAY 1000045 

# BR 
9.81192113E-01 
1.88078869E-02 

DECAYTREE 1000035 

# BR 
1.00000000E+00 

DECAY 1000035 

# BR 
1.00000000E+00 

DECAYTREE 1000025 

# BR 
1.00000000E+00 

DECAY 1000025 

# BR 
1.00000000E+00 

DECAYTREE 1000037 

# BR 
7.89761843E-01 
2.10238157E-01 

DECAY 1000037 

# BR 
7.93795442E-01 
2.06204558E-01 

DECAYTREE 1000024 

# BR 
1.00000000E+00 

DECAY 1000024 

# BR 
1.00000000E+00 



3.87225368E+02 
9.76694000E+02 
1000045 



# 
# 



chil + 
chi2+ 
01 



1.68667313E 
NDA ID1 
2 1000037 
2 1000024 
64835565E-01 # 
NDA ID1 
2 1000037 
2 1000024 

2.33679246E+00 
NDA ID1 
2 1000024 
40601180E+00 # 
NDA ID1 
2 1000024 

2.30613500E+00 
NDA ID1 
2 1000024 
38181900E+00 # 
NDA ID1 
2 1000024 

2.98577444E+00 
NDA ID1 
2 1000023 
2 1000022 
06171556E+00 # 
NDA ID1 
2 1000023 
2 1000022 

2.75859651E-03 
NDA ID1 
2 1000022 
78599108E-03 # 
NDA ID1 
2 1000022 



# S 
ID2 

-24 
-24 

S on NLO 

ID2 

-24 
-24 

# chi04 
ID2 

-24 



# BR(S -> 

# BR(S -> 



# BR(S 

# BR(S 



chi2+ W-) 
chil+ W-) 



chi2+ W-) 
chil + W-) 



#BR(chi04 -> chil+W-) 



chi04 on NLO 
ID2 

-24 #BR(chi04 -> chil+W-) 

# chi03 
ID2 

-24 #BR(chi03 -> chil+W-) 
chi03 on NLO 
ID2 

-24 #BR(chi03 -> chil+W-) 

# chi2+ 
ID2 

# BR(chi2+ -> chi02 W+) 

# BR(chi2+ -> chiOl W+) 
NLO 



24 
24 

chi2+ on 
ID2 

24 
24 



# BR( chi2+ 

# BR( chi2+ 



chi02 W+) 
chiOl W+) 



Block NLOPAR # Renormalization parameters 



# chil+ 
ID2 

24 #BR(chil + 
chil+ on NLO 
ID2 

24 #BR(chil + 



-> chiOl W+) 



-> chiOl W+) 



# UV divergence : 1 = only UV div . 
0.00000000E+00 # UV divergence : Delta 



parts 



1.00000000E-05 
9.11870000E+01 
1.00000000E+05 

1 # Special 
1 # Special 
1 # Special 



# IR divergence : Photon regulator mass 

# Renormalization scale : Q for NLO calc 

# Gauge dependence : Xi 
choice of gauge for photon 
choice of gauge for Z boson 
choice of gauge for W boson 



# NLO masses used for process: LO, 1 NLO masses 
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